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'nJ" ■ Abstract. We model the effects of a large number of zero and near-zero modes 

in the QCD partition function by using sparse chiral matrix models with an 
emphasis on the quenched topological susceptibility in the choice of the measure. 
At finite temperature, the zero modes are not affected by temperature but are 
^\ . allowed to pair into topologically neutral near-zero modes which are gapped at 

0^ ' high temperature. In equilibrium, chiral and U(l) symmetry are simultaneously 

restored for total pairing, evading mean-field arguments. We analyze a number 
of susceptibilities versus the light quark masses. At the transition point the 
topological susceptibility vanishes, and the dependence on the vacuum angle 9 
(j^ " drops out. Our results are briefly contrasted with recent lattice simulations. 



X INTRODUCTION 



The current theoretical resolution of the U(l) problem in QCD relies on 
the assumption that the QCD vacuum supports a finite topological suscepti- 
bility [1-3]. This assumption is supported by current lattice simulations [4], 
anomalous Ward identities [5], chiral effective Lagrangians [6] and canonical 
quantization [7], although there are questions in covariant quantization [8]. 

In this paper, we will adopt the current view and proceed to analyze what 
happens to the U(l) problem at finite temperature. At infinite temperature, 
both the anomaly and topological effects become negligible, so that the U(l) 
symmetry is effectively restored. In this limit there is an exact chiral and 
U(l) degeneracy modulo quark masses. The question then is what happens 
at finite temperature? Does the U(l) restoration coincide or differ from the 
conventional chiral restoration [9]? 
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Recently, a number of lattice simulations [10-13] and model calculations [14] 
have attempted to answer this and other questions with somewhat opposite 

conclusions. It is our purpose in this letter to try to address some of these 
issues in a lattice motivated matrix model, focusing on the interplay between 
zero and near-zero modes, the effects of light quark masses and the impor- 
tance of the thermodynamical limit. In many ways our analysis will parallel 
instanton-like calculations in a solvable context, with interesting lessons for 
these calculations as well as lattice simulations. Indeed, one of the main thrust 
of the present letter is to provide a minimal framework for a model analysis 
of current lattice results. 

In section 2, we motivate the use of a class of chiral matrix models by re- 
viewing some recent lattice calculations. In section 3, we discuss the saddle 
point results following from the present model under some generic assump- 
tions on the interplay between zero and near-zero modes. In section 4, we 
address certain aspects of the chiral and U(l) transitions, including a number 
of susceptibilities. In section 5, we comment on a number of recent lattice 
simulations in light of our results. Our conclusions are in section 6. 
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FIGURE 1. Lattice results from Rcf. [10] for quenched f3 ~ 6.2 (left) and staggered, 

Nf = 2, (3 = 5.55, and niqa = 0.00625 (right) lattices for topological sectors (crosses), 1 
(open circles) and 2 (full boxes). 



Lattice Motivation 

Recently Kogut, Lagae and Sinclair [10] have studied the chirality content 
Tn = (nlj^ln) of the low-lying quark eigenstates A„ for staggered fermions 
with D\n) — Xn\n). Their results after coohng are displayed in Fig. 1 in 
the (r. A) plane for f3 = G/g"^ = 5.55 and 6.2, on a 16^ x 8 lattice. |r| = 1 
in the continuum (about 1/4 on the lattice) corresponds to an eigenvalue 
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with topological charge ±1, while r = corresponds to a non-topological 
eigenvalue. The high temperature configurations (/3 = 5.55) are characterized 

by a depletion in the zero modes (r = 1/4), and an enhancement in the near- 
zero modes (r = 0). Throughout we will refer to the modes with definite 
chirality as zero modes, and those without as near-zero modes. We note that 
at high temperature and in the continuum, the near-zero modes are gapped 
by ±7rr. 



Model 



A simple way to analyze the interplay between the zero and near-zero modes 
around the chiral transition point is through a matrix model. A pertinent ex- 
ample for the fermion matrix was discussed in [15] (and references therein) 
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and the partition function is 



Z[m, 9] = (det D) . 



(2) 



For each fiavor, the entries in the matrices have respectively n, n, n+ and n_ 
elements corresponding to the number of right handed near-zero modes, left 
handed near-zero modes, right handed zero modes and left handed zero modes. 
Here d denotes a diagonal matrix with equal entries, d, the matrices assigned 
to the near-zero modes are square matrices, while the ones assigned to the 
zero modes are rectangular matrices. The fluctuations in the rect angularity 
of the matrices induce the proper U(l) breaking [16]. The hopping between 
zero and near-zero modes is characterized by the overlap matrices T. 

The averaging in (2) is done with respect to the local fluctuations in the 
topological charge x = (^+~^-)) with a Gaussian width flxed by the quenched 
topological susceptibility 



(n. 



2x,V 



(3) 



and a Gaussian measure for the random matrix elements AS, B, with width 
E = 1. The latter is physically tied to the quark condensate |((/tg)| ~ 
(200 MeV)^. For simplicity, we take the width of the Gaussians to be the 
same since the hopping between the low-lying modes may be random enough 
not to distinguish between zero and near-zero modes. The temperature ef- 
fects on the near- zero modes are parameterized by the deterministic and off- 
diagonal entries d. They cause the near-zero modes to be gapped by typically 
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d — ±7rT at high temperature, setting the range of vahdity of the current 
assumptions [15]. The depletion in the number of zero modes caused by an 

increase in the temperature will be discussed below. 

We observe that the columns and rows in the fermion matrix (1) may be 
rearranged to give instead 
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where we grouped the right and left handed modes, respectively. We denote 
the total number of modes (size of the full matrix) by 2A'" = 2n + n+ + n_. We 
note that in the limit mV <tC 1 the matrix model (2) yields sum-rules for the 
quark eigenvalues that are consistent with those discussed in [17], following 
the general arguments in [18] with the identification N = nV, where n is an 
order 1 quantity measuring the density of modes. Physically, the latter sets 
the scale for the fermionic contribution to the energy density. 

Distribution of Eigenmodes 

The matrices D in (1) have (in the chiral limit) a continuous distribution of 
eigenvalues p(A), with a superimposed Dirac delta function at zero virtuality 
A = 0, for fixed topological charge x, 

27Vp(A) = |x|5(A) + 2iVp(A). (5) 

To assess the dependence of Vn on A one can use the fact that A„ and £r„ are the 
real and imaginary parts of the non-hermitian operator D + eil^, in first order 
perturbation theory. Here 15 = diag(l, —1, 1, —1) with each identity assigned 
to its pertinent subspace. Taking e infinitesimally small and rescaling the 
imaginary part of the eigenvalue we obtain the abovementioned dependence^. 

One can now use the chiral structure of D (its anti-commutativity with I5) 
to show that the square 

{D + leUf = - (6) 

is a hermitian operator. This means that the pair A„ + ieVn is cither purely 
real or purely imaginary. It follows that all nonzero eigenvalues (A^ ^ e) 
have vanishing while the |x| topological ones have zero eigenvalue and 
r„ = ±1. In the limit mV > 1 the random matrix model (2) allows for a 
model dependent assessment of the distribution of the low-lying modes of the 
Dirac operator in the infinite volume limit using the methods discussed in [20]. 

It is interesting to note that non-hermitean operators of this form with e = 1 /V^N yield 
a generic distribution of non-hermitean eigenvalues [19] . 
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PARTITION FUNCTION 



Bosonization 



For equal masses, the partition function (2) can be readily bosonized. The 
result for fixed size matrices is 



Here z = me stands for degenerate flavors. The number of near-zero modes 
n, and the number of zero modes n±, fix the size of the matrices in (7). 
However, their distribution is partly fixed by the Gaussian distribution (3), 
while the remaining part is fixed by equilibrium arguments as we now discuss. 



At finite temperature the change in the total number of zero modes 
can be argued generically. Indeed, with increasing temperature the zero 
modes may deplete either by pairing into topologically ncTitral aggregates 
of near-zero modes [21] or screening [22]. For the simplest neutral aggregate 
(molecule [21]), this pairing is reminiscent of the Kosterlitz-Thouless one in 
four-dimensions [23]. 

The chemistry of small neutral aggregates can be described by the proba- 
bility of formation pf and breaking pi,. For molecular arrangements in equi- 
librium, detailed balance implies 



The l.h.s stands for the number of pairs broken. The r.h.s stands for the num- 
ber of pairs formed which is the formation probability, times the probability 
(n_/iV) to find an unpaired negative charge, times the total number of pos- 
itive charges n+. In equilibrium, the number of pairing matches the number 
of breaking. 

We may now calculate the square of {n+ + n_) from the constraint 2N — 
2n + n+ + n_ and subtract x to obtain 



\z+P)(z+P'') + (fY (z+Py+ (z 



n- 



(7) 



Detailed balance 




(8) 



n+n- = {N 



\2 ^2 



(9) 



Therefore n satisfies 




N{l + 6)n 
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X = 



(10) 
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with 5 = Pb/2pf. Since from (3) % ~ \fN we obtain 

„ = iV(l+W5^-^^xVO(/) (11) 

a 

and 

n± = 7V-n±|. (12) 

We are effectively left with a 'fiUing fraction' a and a contribution to the 
topological susceptibility For a ^ 1 we have 5 — (p?, ^ p/). In this 
case, practically all the zero modes arc paired (the unpaired ones arc of order 
l/\fN) with U(l) effectively restored ^. For a — >• we have 5 — >• cxo (p^ 3> p/). 
In this case, all the zero modes are unpaired and U(l) is broken. 



Saddle point analysis 

We insert (11-12) into the partition function (7), and perform a linear shift 
in X, 

X^X-i2N-y (13) 

with the requirement that the term linear in x vanishes in (7). The resulting 
consistency condition (saddle point) reads 

^log|^ + 2iay = (14) 

where 

2 log I I PI2 • 

The parameter y is just proportional to the average topological charge 

{n+-n_)^2Vny, (16) 

while P and in the above equations are the saddle point solutions following 
from the 'action' 

(1-a.) log \z+P\'+a. \og{\z+P\'+d') - iy log - PP^ + — y' (17) 

Z-\-r^ X* 

where the effective 'filling fraction' is 



^) Eq. (11) holds for 1 — a ^ A'' hence in this paper the limit a — > 1 is understood 
always after the thermodynamical limit N ^ oo. 
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Writing out the saddle point equations for P and and subtracting yield 

zP - zP^ = 2iy . (19) 

This suggests the decomposition e~*^P = Q + iy/m, with Q real, being the 
chiral condensate | (g^g) | and satisfying the saddle point equation 

(1-mQ-Q'-^) [{m+Qf+^+d^] = a.d' . (20) 
This equation will be analyzed next. 



RESULTS 

The model is totally specified by (1-3) and (8). The thermodynamical 
limit will be understood as N,V — > oo, with N/V fixed. The parameters 
are: the width of the Gaussian S = 1, the current mass m, the quenched 
topological susceptibility %*, the deterministic entries d, the vacuum angle 9, 
the filling fraction a and the mode-density n — N/V. Generically, the effects 
of temperature cause < d — ttT and < o; < 1. 



Chiral condensate 

The X saddle point equation (14) has a trivial solution, y = for 9 = 0. 
Inserting this back into the equation for the condensate Q and setting m — 
yields 

Q2 ^ 1 (^1-^2 + ^(1-^2)2 + 4(l-«) . (21) 

This result is similar to the one considered by [24], although our physical 
interpretation is different. Indeed, in our case a measures the amount of 
U(l) breaking and follows from the rectangular character of the matrices as 
opposed to the square matrices used in [24]. It is fixed by detailed balance. We 
have ignored the trivial solution with Q = —m, by maximizing the effective 
action (17) 

F/N = -Q' + log{m+Qr + a log . (22) 

For d < 1, the solution with Q ^ sets in independently of a. For a = 1, 
the zero modes pair into near-zero modes, and = 1 — d^ [15,24,25]. This is a 
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U(l) symmetric phase with broken chiral symmetry. A quahtative assessment 
of the range of temperature where this can take place follows by reinstating 
the dimensionful constants, that is c? = ttT < a/S. Hence T < 70 MeV, 
which is outside the range of validity of our model (see above). However, this 
points to the fact that the near-zero modes are sufficiently gapped at already 
moderate temperatures, leaving the zero modes as the only contributors to 
the chirally broken phase. Indeed, at = 1 we have Q = \/l — a, which is 
zero-mode driven. From here on, only the case with d > 1 will be discussed 
unless specified otherwise. 
For a ^ 1, 



Q{m) 



d 



VdF^ 



Vl-a + 0{m) 



(23) 



while for a 



Q{m) = 



m 



d^ 



m 



(24) 



The pairing mechanism suggests an integer 'exponent' S = 1. We recall that 
ior d = a = 1, Q = m}/^ and 5 = 3 which is mean-field [15,24,25]. 



Isotriplet susceptibilities 

A measure of U(l) breaking in the matrix model can be assessed by inves- 
tigating the difference in the tt'^ and aP isotriplet susceptibilities [10,11] 

u = (^qHUr^qqUl^T^qj - (qW^qq^r^q"^ (25) 

and is amenable to the quark eigenvalue distribution p(A) through [12] 

uj^Am'TdX ^^^\ (26) 

For a ^ 1, the matrix model yields Q ^ with a gapped spectrum in the 
chiral limit, hence a; = ^. This observation is similar to the one we made 
in [26] without due care to the U(l) problem as we noted. In general, cu can 
be related to the resolvent 

Specifically, 



^) In fact it would suffice that p(A) vanishes as ~ X°"m'' with a > — 1 and a + b> 1. 
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uj = ^ — - Re G" im) = - m) (28 

m m 



where Q{m) follows from (20). 
For a < 1 we have 



which is to be compared to a; ~ -^1 — a/m for o? = 1. The 1/N corrections to 
(29) are 

+ (...). (30) 

The first term is the contribution of the zero modes in (5) through (26). Since 
X ~ ^fN both contributions in (30) arc subleading in comparison to (29) in the 
thermodynamical limit. These effects may still be present in current lattice 
assessments of uj as we discuss below. 
For a = 1, we have 

uj^^J^y + Oim') (31) 

implying that u flips from 1/m to m? at the transition point. We note that uj ~ 
for (i = 1, a = 1 which is the mean-field result [15]. It is noteworthy 
that only integer 'exponents' are produced by the pairing transition, a point 
in support of some general arguments made in [11]. 



Topological susceptibility 

The topological susceptibility in the matrix model is simply given by 

X..= -^'ogZ = -2|. (32) 

Expanding the consistency equation to linear order in y, we obtain 

y ^ 2a +^ 

^" ^ m(m+Q) 

with a defined in Eq. (15), which gives 

11^1 « , (m+g)2 + rf2 



Xtop X* ~{2mi{mi + Qi) 2(1-q;2) {m+Qy 
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where we have reinstated the flavor dependence. The flrst contribution is the 
quenched susceptibihty, the second contribution is the screening caused by 

the near-zero modes and the unpaired zero modes, and the third contribution 
stems from the paired zero modes. Note that Xtop vanishes not only for mass- 
less quarks but also for maximal pairing with a = 1, as the asymmetry of 
D's become minimal. This happens as Q — > 0, in qualitative agreement with 
recent lattice simulations [4] . 



Pseudoscalar susceptibilities 

The connected and disconnected pseudoscalar susceptibilities associated 
with q^lr-^q may be assessed in a similar way. These susceptibilities were re- 
cently addressed on the lattice [10]. In our case, the disconnected part xf^ 
reads 

Xf = ^ (Trl5-^l5-^) (35) 
N \ im — D im — D / 

and is readily amenable to (34) through 

dis _ 1 / (n+-n_)n _ Xtop .ofiN 
~V\ /~^- ^^^^ 

In the broken phase Xtop is dominated by the the second term in (34) for 
small m, hence xf^ ~ 1/m. As Q ^ 0, the hmits m ^ (chiral) and 
a — >■ 1 (pairing) do not commute. For fixed mass and a ^ 1, xf'"^ ~ (a — 
1) logm/m^ ~ 0, while for a = 1 and m ^ 0, xt'"^ ~ 1/(1 + Q/tu) ~ 1. In 
both cases, xf^ is finite. Note that for d = 1, xf^ ~ m^/'^. 
The connected part Xg""^" follows from the identity [16] 



X.top 




This is the random matrix version of the QCD anomalous Ward identity [5]. 
Hence 

Xf-=(7Vj + l)xf--2^ (38) 

for m > 0. Again, the connected part of the susceptibility is plagued with 
similar ambiguities in the chiral and pairing limits. For a — 1 and m — > 0, 
Xl"""' is finite. 
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Scalar susceptibilities 

The connected and disconnected isosinglet susceptibility associated with q^q 
may be estimated in our case as well, following the lattice conventions [10,11], 

conn ^}_/ rj^ 1 1 \ /ggN 

N\ im-Dim-D/ ^ ' 

and 

drsc ^ 1 / \ \ \ _ J_ / Tr 1 \' . (40) 

N\ im-D im-D/ N^\ im-D/ ^ ' 

Both susceptibihties follow from (20). Specifically, 

= Q\m) = -J—^ (41) 
xr = Q\m) = j^^^. (42) 

for a = 1. This is to be compared with the mean-field result for d = 1, 

■y/^rin _ aj^j Xs'"^ = w?!"^ . The factorized result for the disconnected 

isosinglet susceptibility follows from the absence of correlations in the number 
(n+ + n_). 



Q angle dependence 

In the symmetric phase and for small m, the Q dependence of the free energy 
InZ/y is simple. Indeed, since y ~ m, for a < 1 we may neglect the last term 
in the consistency equation (14) and obtain 

^arctan^ = -e. (43) 



The saddle point equation (20) in the chiral limit can be solved. Defining 
Qi + iy /rrii = \Qi*\e^'^\ the result for each flavor is 

Ql ■ (44) 



m^sin (pi 

where Qj* follows from Q in (21) through the substitution o; — > a* for m = 0. 
Hence 

E'/'^^^' (45) 

i 

mi sin 01 = . . . = mjvy sin 0jv/ ■ (46) 
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These equations are analogous to the zero-temperature equations originally 
derived in QCD [1-3] and more recently in a matrix model [27]. Therefore 

the dependence of the free energy on 6 in the broken phase is the same as the 
vacuum one. The temperature dependence is only implicit through Q*. 

As approaches zero at the critical point and in the chiral limit, the 
dependence on 6 changes. For small m, we may no longer neglect the last 
term in the consistency equation (14) as it diverges. Geometrically the line 
that intersects the curves of the arctan^s becomes nearly vertical so that y 
is for all purposes regardless of the value of the 9 angle. This extends the 
result Xtop = obtained earlier at ^ = to ^ ^ 0. 

The fact that the free energy no longer depends on 9 at the critical point 
and beyond, may be traced to the occurrence of a non-analytic term \x\ in the 
partition function. Indeed from (10) and for a = 1 

^ = ^-^1x1- (47) 

Inserting this into the partition function (7), we obtain 

e^exe-6|xle-2^ (48) 

with b a positive factor stemming from (7). Performing the integral/sum over 
X gives a vanishing contribution to the free energy \ogZ/V. Specifically, the 
sum over x is for a range of parameters well approximated by 

21^^ (49) 

which gives zero contribution to the free energy. This is a direct consequence of 
the total quenching of the topological fluctuations in the paired conflgurations 
of zero modes. The simultaneous restoration of chiral and U(l) symmetry at 
finite temperature yields a symmetric phase that preserves strong CP. 



COMPARISON TO LATTICE 

In a first lattice study by Bernard et al. [11], chiral symmetry restoration 
was found to precede the U(l) restoration. Their analysis relied on gauge 
configurations at fixed lattice spacing a ~ l/6Tc ~ 0.25 fm [11] for Nt = 6. 
Since finite volume effects were not investigated, it may be that the small U(l) 
breaking effects detected in these simulations through a lattice measurement of 
cu for staggered fermions are of the type (30). However, simple estimates based 
on their numbers appear to be on the larger side of their reported results [20]. 
As we already noted, the pairing mechanism supports integer 'exponents' for 
LU, a point sought in [11]. 
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In a second lattice study by Kogut et al. [10], the low-lying quark eigenvalues 
of the staggered Dirac operator where investigated. Their analysis shows that 
the disconnected isosinglet susceptibility x^**, decreases but remains finite 
in the high temperature phase. The finite result was shown to follow from 
the eigenmodes with finite chirality (topological). The conclusion was that 
the U(l) symmetry was not restored in the symmetric phase, although again 
finite volume effects were not investigated. In the present matrix model, we 
have observed that xt''^ remains finite in the chiral and U(l) symmetric phase 
for d > 1, when the thermodynamical limit is carried. Also, we have noted 
an ambiguity in the limits m, a — > 0, 1, suggesting that the cooling procedure 
may be subtle while carrying the chiral limit. Indeed, lattice cooling affects 
the "filling fraction" a. 

In a third lattice study by Chandrasekharan et al. [12], the chiral condensate 
and u) were calculated using also staggered fermions for fixed /3 = 5.3 and 
Nt — A. Although their results were found to be consistent with those of 
Bernard et. al. [11], they concluded that the anomalous effects were small, 
hinting at the possible restoration of U(l) in the symmetric phase. Although 
their conclusions are closer to ours in spirit, they differ in content since their 
small value of uj was obtained from a linear extrapolation in the current quark 
mass, as opposed to a quadratic extrapolation suggested by our results. Also, 
we have observed that the ^-dependence drops in the symmetric phase in 
distinction to a general assumption they made. 

In a fourth lattice study by Vranas et al. [13], lattice simulations with do- 
main wall fermions were carried at Nt — 4. It was found that the high temper- 
ature phase preserves chiral symmetry with a small amount of U(l) breaking, 
although with a somehow heavier pion mass. The method preserves fiavor 
symmetry and incorporates the effect of the anomaly at every stage of the 
simulation. It is indeed encouraging that the results of these simulations are 
the closest to ours. 

CONCLUSIONS 

We have used a simple matrix model to analyze the interplay between zero 
and near-zero modes at finite temperature. While the model finds its mo- 
tivation in the lattice results described above, it was originally argued from 
an NJL model with U(l) breaking [15]. At finite temperature, the pairing 
mechanism at work in the zero mode sector is reminiscent of the one orig- 
inally suggested in the context of instantons [21]. The present model is by 
no means exhaustive as additional effects, e.g. Debye screening, have been 
omitted. Their consideration goes beyond the scope of this work. 

This notwithstanding, our results indicate that chiral and U(l) symmetry 
are simultaneously restored for maximum pairing of zero modes. Although the 
chiral condensate receives contribution from all low-lying modes, its depletion 



to zero requires that the zero modes are paired and the near-zero modes are 
gapped. A simple estimate shows that the near zero modes are substantially 

gapped at moderately low temperatures, suggesting their early decoupling. 
This rules out the possibility of a U(l) restoration prior to a chiral restoration, 
and suggests that both symmetry restorations occur simultaneously. 

The transition by pairing the topological charges is followed by a number 
of observations regarding the topological, scalar and pseudoscalar susceptibil- 
ities for small current quark masses. In particular, integer 'exponents' were 
observed in contrast to the fractional exponents expected from general univer- 
sality arguments. These susceptibilities have been extensively studied on the 
lattice. Our comparison with the most recent lattice simulation using domain 
wall fermions is very encouraging, although some improvements regarding the 
extrapolation to zero qiiark mass and finite volume effects are still warranted 
in the staggered simulations. In many ways, our results should benefit the 
more complex instanton calculations when they become available. 

Finally, we have shown that in the symmetric phase the topological sus- 
ceptibility vanishes in the thermodynamical limit. As a result, the partition 
function develops a non-analyticity in the net topological charge that causes 
the symmetric phase to be CP even whatever the vacuum angle. While ad- 
mittedly this is a result of the present matrix model, it should be interesting 
to see whether it carries to QCD in the infinite volume limit. 
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